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Abstract 

In this paper, we establish a dynamic game to allocate CSR (Corporate Social Responsibility) 
to the members of a supply chain. We propose a model of a three-tier supply chain in a 
decentralized state which includes a supplier, a manufacturer and a retailer. For analyzing 
supply chain performance in decentralized state and the relationships between the members 
of the supply chain, we use a Stackelberg game and consider in this paper a hierarchical 
equilibrium solution for a two-level game. In particular, we formulate a model that crosses 
through multi-periods with the help of a dynamic discrete Stackelberg game. We obtain an 
equilibrium point at which both the prohts of members and the level of CSR taken up by 
supply chains is maximized. 
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1 Introduction 

In recent years, companies and hrms have been showing an ongoing interest in favor of CSR. This 
is mainly because of increasing consumer awareness of several CSR issues, e.g. the environment, 
human rights and safety. In addition, the hrms are also forced to accept CSR due to government 
policies and regulations. Recently CSR has gained recognition and importance as held of research 
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field [HI US]- However, the research field still lacks a consistent definition of CSR and this has 
been the center of discussion since several decades. Dahlsrud [S] presented an overview of different 
dehnitions of CSR and summarized the number of dimensions included in each dehnition. There 
is a positive correlation between CSR and proht [191120] • Moreover, CSR is an effective tool for 
supply chain management, for coordination, purchasing, manufacturing, distribution, and market¬ 
ing functions [TT]. According to previous studies, the long term investment on CSR is beneficial 
for a supply chain. Furthermore, a sustainable supply chain requires consideration of the social 
aspects of the business [23]. Carter et ah [6] established an effective approach and demonstrated 
that environmental purchasing is significantly related to both net income and cost of goods sold. 
Carter et ah [7] also pointed towards the importance of CSR in the supply chain, in particular 
the role played by the purchasing managers in socially responsible activities and the effect of these 
activities on the supply chain. Sethi izn introduced a taxonomy in which a firm’s social activi¬ 
ties include social obligations as well as more voluntary social responsibility. And, Carroll mi 
developed a framework for CSR that consists of economic, legal, and ethical responsibilities. 

The members of a supply chain take their decisions based on maximizing their individual net 
benefits. In addition, when they need to accept CSR; this situation leads to an equilibrium status. 
Game theory is one of the most effective tools to deal with this kind of management problems. 

A growing number of research papers use game theoretical applications in supply chain manage¬ 
ment. Cachon et al. [3] discuss Nash equilibrium in noncooperative cases in a supply chain with 
one supplier and multiple retailers. Hennet et al. [13] presented a paper to evaluate the efficiency 
of different types of contracts between the industrial partners of a supply chain. They applied 
game theory methods for decisional purposes. Tian et al. [2l] presented a system dynamics model 
based on evolutionary game theory for green supply chain management. 

In this paper, we consider a discrete time version of the dynamic differential game. The optimal 
control theory is the standard tool for analyzing the differential game theory [TB]. There are two 
different types of information structures in a differential game, open-loop and feedback information 
structures. In an open-loop strategy, the players choose their decisions at time t, with information 
of the state at time zero. In contrast, in a feedback information structure, the players use their 
knowledge of the current state at time t in order to formulate their decisions at time t HI]- We 
formulate a model and study the behavior for decentralized supply chain networks under CSR 
conditions with one leader and two followers. The Stackelberg game model is recommended and 
applied here to find an equilibrium point at which the profit of the members of the supply chain 
is maximized and the level of CSR is adopted in the supply chain. We develop an open-loop 
Stackelberg game by selecting the supplier as the leader and both the manufacturer and the retailer 
as the followers. Using this approach, the supplier as a leader, can know the optimal reaction of 
his followers, and utilizes such processes to maximize his own profit. The manufacturer and the 
retailer as followers, try to maximize their profits by considering all the conditions. Our model 
has two levels, at the first level the manufacturer is the leader and the retailer is the follower and 
we hud the equilibrium point. At the second level, we consider the supplier as the leader and the 
manufacturer as the follower. In fact, we substitute the response functions of the follower into the 
objective function of the leader and we find the final equilibrium point. We propose a Hamiltonian 
matrix to solve the optimal control problem to obtain the equilibrium in this game. The paper is 
organized as follows: Section 2 is devoted to the problem description and assumptions. Objective 
functions, constraints and solution of the game are illustrated in Section 3. A numerical example 
is shown in Section 4 and we close with a conclusion in Section 5. 
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2 Problem Description and Assumptions 


We consider a three stage Stackelberg differential game which has three players playing the game 
over a fixed hnite horizon model. This model is a three-tier, decentralized vertical control supply 
chain network (Figure 1). All retailers and suppliers at the same level make the same decision. 
Therefore, consequently the model has only one supplier, one manufacture and one retailer. The 
simplihed model is shown in Figure 2. 

The dynamic game goes through multi-periods as a repeated game with complete information. 
This model has a state variable and control variables like any dynamic game. We dehne the state 
variable as the level of social responsibility taken up by companies, and the control variables are 
the capital amounts invested while fulfilling the social responsibility. Specihcally, all of the social 
responsibilities taken up by the firm j at period t can be expressed as the investment /^. We 
suppose that Xt evolves according to the following rule: Xt+i = f{It,Xt). 



Figure 1: Three-tier supply chain network 


More specihcally we have the following assumptions: 

The function Bt{xt) = 5xt represents social beneht which is proportional to social responsibility 
taken up by the supply chain system [T]. 

The function T* = r/t[l -|- 9{It)\ measures the value of the tax return to the members of the 
supply chain [lO]. Both r and 6 are tax return policy parameters. Specihcally, r is the rate of 
individual post tax return on investment (ROI), and 6 is the rate of supply chain’s post tax return 
on investment (ROI). 

The market inverse demand is P^{qt) = a — bqt [T7] . 


3 







































QJ 


0 0 0 

f —^ 

-A 



M 








C —’ 






5: Supplier 








R\ Retailer 



( R 


( Period t ] I Period t+1 I 

^ I 


Figure 2: Simplified model of three-tier supply 
chain network 

The accumulation of the level of social responsibility taken up by the firms is given by Xt+i = 
axt + ■ 

Here, /3i is the rate of converting the supplier’s capital investment in CSR to the amount of CSR 
taken up by the supply chain; A is the rate of converting the manufacturer’s capital investment 
in CSR to the amount of CSR taken up by the supply chain and A is the rate of converting the 
retailer’s capital investment in CSR to the amount of CSR taken up by the supply chain as well 

m- 


2.1 The General Framework 


He et ah [12] illustrate an open-loop Stackelberg differential game model over a fixed hnite horizon 
time as detailed in the following: 


The follower’s optimal control problem is: 


Max 


r(-) 


subject to the state equation 




X{t)=F{X{t),w{t),r{t)), 
X(0) = Xo. 


( 2 . 1 ) 


( 2 . 2 ) 


where the function F represents the rate of sales, p is the followers’s discount rate, and Xq, is the 
initial condition. The follower’s Hamiltonian is 


//h(V, r, Ah, w) ^ 7Th(A:, w, r) + \hF{X, w, r), 


(2.3) 


where Ar is the vector of the shadow prices associated with the state variable X; and it satishes 
the adjoint equation 


Ar — pAr — 


dHR{X,r, Ar,w) 
dX 


Ar{T) = 0. 


(2.4) 
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The necessary optimality condition for the follower’s problem satishes 


dH 


R 


= 0 


dr dr 

We derive the follower’s best response r*{X,w, Xr). 
The leader’s problem is 


d7rR{X,w,r) dF{X,w,r) 

-H - — U. 


dr 


(2.5) 


Max^(.)| = / e '^^7iM{x,w,r{x,w, 


X = F{X,w,r{x,w,XR)), 

Xo{0) = Wo, 

dHR{x, r{x, w, Xr), Xr, w) 


( 2 . 6 ) 


^R — P^R — 


dx 


, ^RiT) — 0, 


where /r is the leader’s disconnt rate and the above differential equations are obtained by substi¬ 
tuting the follower’s best response r*{X,w, Xr) into the state equation and the adjoint equation 
of the follower, respectively. We formulate the leader’s Hamiltonian as follows 


dHR{X,r{X,w,XR),XR,w) 
dX 


Hm = Xr, w, r{X, w, Xr), Xm, ^) + XF{X, w, r{X, w, Xr)) - fi 

( 1 A. 

(2.7) 

where Xm and /r are the shadow associated with X and Xr, respectively, and they satisfy the 
adjoint equations 


Xm — pX 


M 


— pXm 


dHMjX, Xr, w, r{X, w, Xr),Xm, w) 
dX 

d'KM{X,w,r{X,w,XR)) 

dX 


^ dF{X,w,r{X,w,XR)) d‘^HR{x,r{X,w, Xr), Xr,w) 
Am-- P - 


ip = pcp- 


dX ^ 

dHM{X, Xr, w, r{X, w, Xr), Xm, p) 


( 2 . 8 ) 


dX 


R 


dF{x,w,r{X,w,XR)) d‘^HR{X,r{X,w, Xr), Xr,w) 
= pip — Xm - Tp: - P - 


dX 


R 


dXdX 


where Xm(T) = 0 and (^(0) = 0 are the boundary conditions. 

We apply the algorithm of the above general model as part of our model. 


2.2 Notations and Definitions 

To facilitate the model, certain parameters and decision variables are used. 


Table [U shows notations and dehnitions that we use in our model. 
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Variables 

t 

T 

Qt 


b 

Xt 

rM 

rR 

jM 

jR 

B^{xt) 

T^xt) 

T^ixt) 

T^ixt) 

jM 

If 

d 

d 


5 

? 

? 

A 


e 

A 


A 


A 


Period t 

Planning horizon 
Demand quantity at period t 
Market potential 
Price sensitivity 

State variable, degree of taking SR 
Hamiltonian function of the supplier 
Hamiltonian function of the manufacturer 
Hamiltonian function of the retailer 
Objective function of the supplier 
Objective function of the manufacturer 
Objective function of the retailer 
Social beneht of the supplier 
Social beneht of the manufacturer 
Social beneht of the retailer 
Tax return of the supplier 
Tax return of the manufacturer 
Tax return of the retailer 

The amount of investment done by the manufacturer 

The amount of investment done by the supplier 

The amount of investment done by the retailer 

The percentage of investment of the supplier payoh 

The percentage of investment of the manufacturer payoh 

Price of the supplier’s raw material 

Price of sold product by the retailer 

Parameter of the supplier’s social beneht 

Parameter of the manufacturer’s social beneht 

Parameter of the retailer’s social beneht 
Quantity discount parameter of the price of raw material 
Deteriorating rate of the level of current social responsibility 
The rate of individual post tax return on investment (ROI) 

The rate of supply chain’s post tax return on investment (ROI) 

The rate of converting the supplier’s capital investment in CSR 

to the amount of CSR taken up by the supply chain 

The rate of converting the manufacturer’s capital investment in CSR 

to the amount of CSR taken up by the supply chain 

The rate of converting the retailer’s capital investment in CSR 

to the amount of CSR taken up by the supply chain 

Table 1: Notations and Dehnitions. 


3 Objective Functions and Constraints 


The objective functions are made to depend on the control vectors and the static variable. The 
members of the supply chain attempt to optimize their net profits, which includes minimizing the 
cost of raw materials and investment in social responsibility, and maximizing sale revenues and 


6 




benefits from taking social responsibility as well as tax returns. Thus, the objective function of 
the supplier is 


= E - <=''< + -Bf(i,) + r/(/f, /,) - /f + dif< 

t=l 

= '£wg,- cq, + fa? + r/f [1 + 9(/f + /" + /«)] - if + diff 

t=l 

where, is the price of the supplier’s raw material. P/ = w. Bf[xt) is the social beneht of the 
supplier, 5 is the parameter of the supplier’s social beneht and is the tax return of the 

supplier. Similarly, the objective function of the manufacturer is 

jM ^ Pf\q,)q, - Pfq, + I,) - /" + dif 

t=l 

T 

= ~ — wqt + + rl ^(1 + 0(jf ++ 1^)) — + dl ^, 

t=i 

where P^{qt) is the retail price of the product of the manufacturer. B^(xt) is the social beneht 
of the manufacturer, 6 is the parameter of the manufacturer’s social beneht and , It) is the 

tax return of the manufacturer. 

The objective function of the retailer is 

jR = J2 PtSt - Pt^{qt)qt + B^{xt) + h) - if 

t=i 

T ^ 

= Y^zq,-(a- bqfq, + fa? + r/«(l + 9(/f + if + if)) - if, 

t=l 

where P^ is the price at which the retailer sells the product to the consumer. P^ = Z. B^{xt) is 

the social beneht of the retailer, 6 is the parameter of the retailer’s social beneht and It) 

is the tax return of the retailer. 


3.1 Mathematical Model: Level One 

At this level, we establish a Stackelberg game between the manufacturer as the leader and the 
retailer as the follower. To calculate the equilibrium at this level, hrst we calculate the best 
reaction function of the retailer, then we determine the manufacturer’s optimal decisions based on 
the retailer’s best reactions. 

Since our dynamic diherential game is an optimal control problem, we can apply the Hamiltonian 
function to hnd the equilibrium of the game IZD. 

Suppose the time interval is [1,T]. For any hxed and I^ the retailer solves 

T 

arg max P,% - Pt^{qt)qt + B^{xt) + T^I^, h) - I^, 

F t=i 
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subject to Xt+i = axt + (5ilf + (^ 2 !^ + 

We define the retailer’s Hamiltonian for fixed and as 

By using the conditions for a maximization of this Hamiltonian, we compute: 


„ 1 - - r«(/" + If) - T 




2«t 


(3.1) 


The equation of which depends on and , says that for any given strategy of Jf and 
there is a unique optimal response 


R 


Xt+l = 




— axt + Pilf + ^2^^ + Psif' 5 


and by substituting (13.Ih in (13.21) . we obtain 


Xt+i — (A — 1^3/2)/^ + (/52 — ^3/2)1^ + ttXt + 


M , , « 1 - ^t+i/^3 - r 


2x6 


We also have 


P = ^ = 2?^:, + aPi,. 

The above sets of equations define the reaction function of the retailer. 

For any fixed if the manufacturer solves 

arg max ^ {qt)qt - Pfqt + B^{xt) + h) - I^ + dif, 

jM f ^ 

‘ t=i 

subject to Xt+i = axt + jdilf + + jd^lf- 


(3.2) 


(3.3) 


(3.4) 


Now, we substitute the value of If in fl3.1l) into and we obtain 


= (a - bqt)qt - wqt + 6xf + 


t 9 — d 


/f If + 


r - 1 - - d 


- jM _ 

' * 26' 


(3.5) 


The Hamiltonian function of the manufacturer for fixed if is 


Hf^ = jf^ + Pfffxi+,)+ufPf), 


(3.6) 


consequently, we can obtain the unique optimal response of the follower from the equations as 
follows. 













dH, 


M 


dL^ 


= r(l + 0((l-r0)jf + J/"))-l + 


1 - PsPt+l - T -Td 


+ (/^2 — /^ 3 / 2 ) Pt 


M 

t+1) 


and we get 


tM ~i'^It + PsPt+i + (1 + d — r) 

- P3/2)Pif, 

(3.7) 

* 2x9 

t6 

Other constraints are 



P)ljM 

Pt^ = ^ = 

P 2dut, 

(3.8) 

^ dp, 

Ht+I ar>R yi! o o 

pftiPi 

O z, +«“i- 

(3.9) 


t+i 


The equation of depends on , and we obtain the hnal equilibrium in the next section, at 
level two. 


3.2 Mathematical Model: Level Two 

At the previous level, the manufacturer’s optimal function was calculated by using a reaction 
function of the retailer. At this level, the reaction functions of two followers (retailer and manu¬ 
facturer) are inserted into the objective function of the leader (supplier), and we can hnd the hnal 
equilibrium point. 

The problem facing the supplier is simply given by 


T 


argmax V Pfqt - cqt + Bf{xt) + Tf (/f, h) - if + dlf^, 
It 


t=i 


subject to xt+i = axt + (Iilf -f ^ 2 ^*^ + Id 2 .lf- 


The Hamiltonian function of the supplier is dehned by 


Hf = y + ). (3.10) 

Substitute (13.ip and (13.71) into (I3.10p . we get the value of Jf, Xt+i and Ht+i 



(3.11) 


9 

















therefore 


s _ (/^3/2 + [^ 2 - 2/3i) s 




tO 


pb 

t+1 


{1^2 - /53/2) j^M , /^3 pi? 

t+1 “T 


t6 


pM , 

t+1 


2r0' 


3 — 3r — c/ + 2(i 
2^0 ■ 


(3.12) 


We have 


therefore we obtain 


Xt+l = 


dHf 

’ 


(3.13) 


Xt+l — OiXt + (/?! — / 32/2 — (3^/4,)!^ + 

(2/32 ~ ( 33 )(1 ~ X + (3) + 2 / 33(1 — t) 


2/32 + / 33)(/32 -/33 / 2 ) 


2r0 


dM 

t+1 


, (2/32/33) - (3/3|) 

“T , ^ -i V 


4r6' 


t+i 


(3.14) 


+ 


At6 


We also have 


/^i+l 


And we obtain 


_ {(32 - 133/2) ^ (/32 -/33 /2)(/33 - 2/32) 

^ - ai". ^ 

(fe - ft/2)rf 

tO 


(3.15) 


Pf = ^ = 25xt + «p4i + 2(5n4. (3.16) 

Since we use open-loop information, the structure variables depend on the time variable and the 
initial state variables. The xi is given initial parameter, ui = 0 and pi = 0. Furthermore, the 
boundary condition are P^i= 0, P//li=Q and Pi^i= 0. 


3.3 Augmented Discrete Hamiltonian Matrix 

In this section for solving the optimal control problem formulated in Section 3.1 and 3.2, we 
chose an algorithm given by Medanic and Radojevic which is based on an augmented discrete 
Hamiltonian matrix [18]. First, we assume 


Xt+l 


■ A 

B ' 

Xt 

+ 

■ D ■ 


Axt P BPt+i P D 

Pt 


C 

D 

. Pt+i . 

E 


Cxt P APt_|_i -|- E 



Xt+l 


' Pf+i ' 

where Xt+i = 

Ut+l 

and Pt+i = 

- 1 
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A, B, and C are 3x3 matrices, and D and E are 3x1 matrices, such that 



Xt+l 




Oil 

012 

Ol3 




Xt+l = 

hf+1 

= 

Axt + i?Pt+i + D 

= 

^21 

O 22 

O 23 


hi 

+ 


Ml+1 




^31 

O 32 

O 33 


Ml 



bn 

^21 

b 


'31 


012 

^22 

^32 


0l3 

bl3 

^33 


di 

d2 

d3 


o-nXt + CLufJ't 

T u.23«'l T 1^21-' t+1 “T '^22-' t+l T 1^23-1 t+1 T u-2 
0‘3lXt + fl32/it + a33'Ol + fesi-Pt+l + ^32-fi+i + &33-ft+l + ^3 _ 


-T ai3Ut -f- unJ^t+i -I- 1*12-'t+1 -r Ci3-r^_,_i -t- ai 
0-2lXt + 0,22fJ‘t + 0,23Ut + 621-Pt+l + &22-ft+l + &23-ft+l + C?2 

, „ 7 -i,S , 7 ->M , 7 _ 7 ^,S , „7 



r ps' 1 

^t+1 


jdM 

^t+1 


PSr^ J 


(3.17) 



■ 1 ■ 


■ 0 ■ 

The boundary conditions are Xi = 

0 

and Pt +1 = 

0 


0 


0 



Oil 

O 12 

O 13 

A = 

021 

O 22 

O 23 

where 

031 

O 32 

O 33 

Oil = 

a. 

O 12 = 

0 , 


■ &11 

^12 

bi3 

B = 

^21 

^22 

^23 

where 

^31 

^32 

^33 


bn = 


(/3i — /32/2 — /33/4)(—2/3i + /32 + /33/2) 

~e ' 


(A - /52/2 -/33/4)(/32 -/53/2) (/32 - /33/2)(-2/32 +/^s) 

O 12 — -;;-r 


tO 


‘ItB 


, (A-/92/2-/33/4)(/33) , (2/32/33)(-3/3|) 

-2^0-+-4^0-• 

{1^2 — /33/2)(2/3i — 3/^2 + /?3/2) 


&21 — 


^22 — 
^23 = 


2t9 

W 2 - 

2 r 0 

Id3{(d2 — /53/2) 


&31 — 
^32 = 


4r0 

/^3(~2/3i + /32 + /33/2) 
4^0 ■ 

'/^l + 3/2/33 (/32 — /33/2) 
2 r 0 
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, where 


hs — 



D = 


di 

d2 

dz 


di = 


(/3i-/32/2-/33/4)(3-3r-c? + 2c?) (2/32 

2 r 0 ^ 


/^3)(1 — T + d) + 2 / 33(1 

At 6 



(^2 


{-(32 - l3z/2){Qd - h - 3r + 3) 

4^0 


— (3z{ — 7d + 2(i — T + 1) 

- 

^ 8t9 

Similarly, we can get the values of the matrices C and E 



' pf ■ 


Cll 

C12 

cl3 ■ 


Xt 


Oil 

0-12 

013 


r p>5 1 

m+i 


ei 

Pt = 

tdM 

= 

C2I 

C22 

c23 


Pt 

+ 

0-21 

O 22 

O 23 


jdM 

^t+1 

+ 

62 


. -P" . 


. ^31 

C32 

c33 


. . 


O 3 I 

O 32 

O 33 


P5i J 


. ®3 . 


CuXt + Ci2/it + Ci3Mi + + auP^li + Q-izPt^i + Ci 

C2lXt + C22fJ't + C2zUt + a21-Pi+i + a22-ft^i + a23-Pf+l + ^2 

cziXt + C 32 /it + czzUt + aziPf_^i + az 2 Pi+i + 033-^4+1 + cs 



2(5 2? 25 


■ 0 ■ 

Therefore C = 

25 0 25 

and E = 

0 




0 


25 0 0 




3.4 Resolution 

The above problem is solved by the sweep method [2], by assuming a linear relation between p/ 
and Xt 


Pk ^kXk 9 k- 


(3.19) 


Thus, we can compute 


Xk+i = {l 2*2 - BSk+i) ^{Axt - Bgk+i + P>)- 


(3.20) 
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Then by substituting fl3.19p and fl3.20p into the dehnition of Pk+i as given by the augmented 
Hamiltonian matrix, and equating both sides we dually get the difference equations: 


Sk 


C + ASk+i{l2*2 — BSk+i) 


(3.21) 


Qk — ASk+i{l2*2 — BSk+i) ^ + Bgk+1 — D + Agk+i — E. (3.22) 



Xi 


■ 0 ' 


The boundary conditions are xi = 

0 

and Rt +1 = 

0 

. And then St+i = 63*3 and 


0 


0 



Qt+i — Os*!. 


From the boundary conditions we get St = C and gx = E. Once we get the different values of Sk 
and gk by the backward loop, then the values of Xt and pt are computed by a forward loop. And, 
consequently we get the values of Xt, , I^, pf, p^, pf, for all points in time. 


4 Numerical Example 


In this section we provide a numerical example. We run the following numerical simulations with 
mathematica 8. The results presented here are obtained for the following values of the parameters: 
a = 6, w = 3.8, c = 2.4, qt = 100000, d = 0.6, c = 0.00001, d = 0.4, d = 0.4, z = Q, 9 = 0.01, 
and r = 0.2. We set /3i = 0.3, (32 = 0.5 and (3^ = 0.8; and /9i = 0.3, (32 = 0.5. Bt{xt) = 6xt, the 
potential benehts hrms obtain from taking social responsibility, such as increased demand, better 

reputation and so on. We set 6 = 0.2, 6 = 0.2 and 6 = 0.2. We assume that the time horizon is 
T=10. The initial level of social responsibility is supposed to be a;i=l. We draw the results of the 
equilibrium from our model, a three-stage Stackelberg dynamic game. 

The hgure [3] shows the trend of prohts from periods one to ten in a Stackelberg game. JS is the 
supplier’s proht, JM is manufacturer’s proht and JR is retailer’s proht. We compare the prohts 
of the supplier, manufacturer and retailer over a time horizon, hrst while playing the game and 
then, without playing the game. Figure 0] shows the difference in supplier’s prohts when playing 
the game and without playing. JSO is supplier’s proht without playing the game; JS is supplier’s 
proht when playing the game. As in the hrst graph, the second and third one (hgure [5l [ 6 ]) shows 
the diherence in manufacturer’s proht and retailer’s prohts when playing the game and without 
playing. JMO is manufacturer’s proht without playing the game; JM is manufacturer’s proht 
when playing the game. JRO is retailer’s proht without playing the game; JR is retailer’s proht 
when playing the game. Obviously, all of players gain extra proht from playing the games. Figure 
[7] compares the cumulated prohts of the member’s of supply chain, playing game one and without 
playing game. 

In sum, the supplier, manufacturer and retailer are motivated to play the game because their 
respective beneht increases and the supplier as the leader in the game earns more beneht than the 
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followers. Of course, this result is obtained with a very specihc dynamic game model. Another 
one may give different results. 



Figure 3: Prohts of supplier, manufacturer and 
retailer. 
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Figure 4: Comparison of the supplier’s proht, 
playing game one and without playing any game. 
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Figure 5: Comparison of the manufacturer’s 
proht, playing game one and without playing any 
game. 
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Figure 6: Comparison of the retailer’s profit, 
playing game one and without playing any game. 
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Figure 7: Comparison of the cumulated profits of 
the member’s of supply chain, playing game one 
and without playing any game. 

5 Conclusion 


In this paper we investigated a decentralized three-tire supply chain consisting of supplier, manu¬ 
facturer and retailer with the aim of allocating CSR to members of the supply chain system over 
time. We considered two-level Stackelberg game consisting of two followers and one leader. The 
members of a supply chain play games with each other to maximize their own profits; thus, the 
model used was a long-term co-investment game model. The equilibrium point in a time horizon 
was determined at where the profit of supply chain’s members was maximized and CSR was im¬ 
plemented among members of the supply chain. We applied control theory and used an algorithm 
(augmented discrete Hamiltonian matrix) to obtain an optimal solution for the dynamic game 
model. We presented a numerical example and we found that, the benehts of the player increased 
when they played the game. 
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